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Abstract. In this article, we are concerned with long-time behaviour of solutions to 
a semi-classical Schrodinger-type equation on the torus. We consider time scales which 
go to infinity when the semi-classical parameter goes to zero and we associate with each 
time-scale the set of semi-classical measures associated with all possible choices of initial 
data. We emphasize the existence of a threshold : for time-scales below this threshold, 
the set of semi-classical measures contains measures which are singular with respect to 
Lebesgue measure in the "position" variable, while at (and beyond) the threshold, all the 
semi-classical measures are absolutely continuous in the "position" variable. 



1. Introduction 

1.1. The Schrodinger equation in the large time and high frequency regime. 

This article is concerned with the dynamics of the linear equation 

f ihdtif> h {t,x) = H{hD x )if> h (t,x), {t,x)eRxT d , 

on the torus T d := (M/27rZ) d , with H a smooth^] real- valued function on (M. d )* (the dual of 
R d ), and h > 0. In other words, if is a function on the cotangent bundle T*T d = T d x (R d )* 
that does not depend on the first variable, and thus gives rise to a completely integrable 
Hamiltonian flow. We are interested in the simultaneous limits h — > + (high frequency 
limit) and t — > +oo (large time evolution). Our results give a description of the limits of 
sequences of "position densities" (t/j,x)| 2 at times th that tend to infinity as h — > + . 
To be more specific, let us denote by the propagator associated with H(hD x ): 

Fix a time scale, that is, a function 

r : R* + — y R* 
h i — y r h , 
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partially supported by grants MTM2010-16467 (MEC), ERC Starting Grant 277778. 

1 For the sake of simplicity, we shall assume that H e C°° (R d ). However the smoothness assumption 
on H can be relaxed to C k , where k large enough, in most results of this article. 
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such that liminf/ l ^ + r h > (actually, we shall be mainly concerned in functions that go 
to +00 as h — > + ). Consider a sequence of initial conditions (uh), normalised in L 2 (T d ): 
||w/ l || i 2^ T d~) = 1 for h > 0, and /i-oscillating in the terminology of [T5| ITT] , i.e.: 

(2) limsup \\l[o,R) (-h 2 A) u h \\ L2 ^ 



ft,— ^0+ 



where 1[o,r] is the characteristic function of the interval [0, R}. Our main object of interest 
is the density IS^Uhf, and we introduce the probability measures on T d 



Uhit.dx) := \S t h Uh(x)\ 2 dx; 

the unitary character of implies that G C (M; V (T d ) ) B To study the long-time 
behaviour of the dynamics, we rescale time by and look at the time-scaled probability 
densities: 

(3) v h (r h t, dx) . 

When t ^ is fixed and Th grows too rapidly, it is a notoriously difficult problem to obtain 
a description of the limit points (in the weak-* topology) of these probability measures as 
h — > + , for rich enough families of initial data Uh- See for instance |31[ I3"U] in the case 
where the underlying classical dynamics is chaotic, the Uh are a family of lagrangian states, 
and Th = h~ 2+e . In completely integrable situations, such as the one we consider here, the 
problem is of a different nature, but rapidly leads to intricate number theoretical issues 
[251 [2H [26]. 

We soften the problem by considering the family of probability measures ([3]) as elements 
of L°° (R; V (T d )). Our goal will be to give a precise description of the set M. (r) of their 
accumulation points in the weak-* topology for L°° (M; V (T d )), obtained as [uh] varies 
among all possible sequences of initial data /i-oscillating and normalised in L 2 (T d ) . 

The compactness of T d ensures that Ai (r) is non-empty. Having v G M. (r) is equivalent 
to the existence of a sequence (h n ) going to and of a normalised, /i n -oscillating sequence 
(u hn ) in L 2 (T d ) such that: 

(4) lim — f [ xi^lSj u^ix)] 2 dxdt = [ [ \ (?) v (t, dx) dt, 
n->+oo r hn J ThnCt J T d J a Jjd 

for every real numbers a < b and every x G C (T d ) . In other words, we are averaging the 

densities |S^u/i(a;)| 2 over time intervals of size Th- This averaging, as we shall see, makes 
the study more tractable. 

If case occurs, we shall say that v is obtained through the sequence (uh n ). To simplify 
the notation, when no confusion can arise, we shall simply write that h — > + to mean 
that we are considering a (discrete) sequence h n going to + , and we shall denote by (uh) 
(instead of (uh n )) the corresponding family of functions. 



In what follows, V (X) stands for the set of Radon probability measures on a Polish space X. 
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Remark 1.1. When the function r is bounded, the convergence of (r^t, •) to an accu- 
mulation point v (t, •) is locally uniform in t. Moreover, v can be completely described in 
terms of semiclassical defect measures of the corresponding sequence of initial data (uh), 
transported by the classical Hamiltonian flow cf) s : T*T d — > T*T d generated by H , which 
in this case is completely integrable. Explicitly, 

<f> a (x,0 := (x + sdH(0,0- 

This is nothing but a formulation of Egorov's theorem (see, for instance, (Hjj. Consider 
for instance, the case where the initial data Uh are coherent states : fix p G (R d ) 
with \\p\\ L 2(^d\ — I? fi x (^O)^o) G IR d x M. d , and let Uh {x) be the 27rZ d '-periodization of the 
following coherent state: 

1 fx-x \ & 
Then (t, •) converges, for every t G R, to: 

b~x +tdH(£ ) (x) . 

When the time scale is unbounded, the t-dependence of elements v G Ai (r) is not 

described by such a simple propagation law. From now on we shall only consider the case 

where — > +oo. 
h — >o 

The problem of describing the elements in Ai (r) for some time scale (r/J is related 
to several aspects of the dynamics of the flow such as dispersive effects and unique 
continuation. In (H |22] the reader will find a description of these issues in the case where 
the propagator is replaced by the semiclassical Schrodinger flow e lhtA corresponding to 
the Laplacian on an arbitrary compact Riemannian manifold (corresponding to H(hD x ) = 
—h 2 A in the case of flat tori). In that setting, the time scale = 1/h appears in a natural 
way, since it transforms the semiclassical propagator into the non-scaled flow e thThtA = 
e ttA . The possible accumulation points of sequences of probability densities of the form 
\e ltA Uh\ 2 depend on the nature of the dynamics of the geodesic flow in the manifold under 
consideration. Even in the case that the geodesic flow is completely integrable, different 
type of concentration phenomena may occur, depending on fine geometrical issues (compare 
the situation in Zoll manifolds |20| and on flat tori [21] [3]). When the geodesic flow has the 
Anosov property, the results in [5] rule out concentration on sets of small dimensions, by 
proving lower bounds on the Kolmogorov- Sinai entropy of semiclassical defect measures. 

1.2. Semiclassical defect measures. Our results are more naturally described in terms 
of Wigner distributions and semiclassical measures (these are the semiclassical version 
of the microlocal defect measures |16[ 152"] . and have also been called microlocal lifts in 
the recent literature about quantum unique ergodicity, see for instance the celebrated 
paper [19]). The Wigner distribution associated to Uh (at scale h) is a distribution on the 
cotangent bundle T*T d , defined by 

(5) / a(x, de) = (u h , Vh (a)u h ) L2m , for all a G C™(T*T d ) } 

JT*1 d 
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where Op h (a) is the operator on L 2 (T d ) associated to a by the Weyl quantization. The 
reader not familiar with these objects can consult the appendix of this article. For the 
moment, just recall that if x is a smooth function on T*T d = T d x (R d )* that depends only 
on the first coordinate, then 

(6) / x(x)i»u h (dx,d£) = / x(x)\u h (x)\ 2 dx. 

JT*T d J1 d 

The main object of our study will be the (time-scaled) Wigner distributions corresponding 
to solutions to (JTJ): 

Whit, ■) := w h Th t . 

o h u h 

The map t i — > Whit,-) belongs to L°°(R; V (T*T d )), and is uniformly bounded in that 
space as h — > + whenever (v.h) is normalised in L 2 fT d ). Thus, one can extract subse- 
quences that converge in the weak-* topology on L°°(R; V (T*T d )). In other words, after 
possibly extracting a subsequence, we have 

ip(t)a(x,£)wh(t,dx,d£)dt — > / / tf{t)a{x,^)ji{t,dx,d^)dt 

! T *jd h >0 J R J T *jd 

for all <p G L^R) and a e C~(T*T d ), and the limit y, belongs to L°° (R; M+ (T*T d ))| 

The set of limit points thus obtained, as (uh) varies among normalised sequences, will 
be denoted by Ai (r). We shall refer to its elements as (time- dependent) semiclassical 
measures. 

Moreover, if (u h ) is /i-oscillating, it follows that y e L°° (R; V (T*T d )) and identity © 
is also verified in the limit : if v (t, •) is the image of \x (t, •) under the projection map 
(x,£) i — > x, then 

b r r-b r 

/ x ( x ) \&h u h ( x ) \ 2 dxdt — > / / x ( x ) A 4 (t, dx, d£) dt, 

a JJ d ' h >0 Ja JT*T d 

for every a < b and every x £ C°° {^ d ) ■ Therefore, Ai (r) coincides with the set of projec- 
tions onto x of semiclassical measures in Ai (r) corresponding to h -oscillating sequences 
(see Pan?]). 

It is also shown in the appendix that the elements of Ai (r) are measures that are 
invariant by the Hamiltonian flow <f) s : 

/ so <j> s (x,g)fi(t,dx,d€) = / a(x,£)n(t,dx,dg), V/x e M (r) ,Vs e R, a.e.t e R d . 

JT*T d JT*T d 

1.3. Regularity of semiclassical measures. The main results in this article are aimed 
at obtaining a precise description of the elements in Ai (r) (and, as a consequence, of those 
of Ai (t)). We first present a regularity result which emphasises the critical character of 
the time scale t/j = 1/h in situations in which the Hessian of H is non-degenerate, definite 
(positive or negative). 



M.+ (X) denotes the set of positive Radon measures on a Polish space X. 
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Theorem 1.2. (1) Ifr^ then M. (r) contains elements that are singular with respect 

to the Lebesgue measure dtdx. Besides, Ai (r) contains all uniform orbit measures of S . 

(2) Suppose Th ~ 1/h or Th ^> 1/h. Assume that the Hessian d 2 H(^) is definite for all£. 
Then 

M (r) C L°° (R; L 1 (T d )) , 
i.e. the elements of M. (r) are absolutely continuous with respect to dtdx. 

Remark 1.3. Theorem \1.2\( 2) applies in particular when the data (uh) are eigenfunctions 
of H(hD x ) , and shows (assuming the Hessian of H is definite) that the weak limits of the 
probability measures \uh(x)\ 2 dx are absolutely continuous. 

Remark 1.4. The conclusion of Theorem \1.2V 2) may fail if the condition on the Hessian 
of H is not satisfied. We give here two counter-examples. 

Fix ui G M. d and take H (£) = £ • to. Let /x be an accumulation point in V (T*T d ) 
of the Wigner distributions (toj J defined in (JSJ) , associated to the initial data (u^) . Let 
/it G M. (r) be the limit of w h Tht in L°°(W;V (T*T d )). Then an application of Egorov's 

theorem (actually, a particularly simple adaptation of the proof of Theorem 4 in [20] ) gives 
the relation, valid for any time scale (t/j) : 

a (x, f ) fj, (t, dx, d£) = / (a) (x, f ) /x (dx, d£) , 

T*T d JT*T d 

for any a G (T*T d ) and a.e. t E R. Here (a) stands for the average of a along the 
Hamiltonian flow <p s , that is in our case 

1 f T , 

(a)(x,£)= lim — / a (x + scj, £) ds. 

Hence, as soon as u is resonant (in the sense of §2. 1\) and [1q = 5 Xo £g> 5^ for some 
{ x o,£,o) £ T*T d , the measure /i will be singular with respect to dtdx. 

It is also easy to provide counter-examples where the Hessian of H is non-degenerate, 
but not definite. On the two-dimensional torus T 2 , consider for instance H(£) = £f — £f ; 
where £ = (^1,^2)- Take for (uh(xi,X2)) the periodization of 

1 ( X\ — X2 

rP 



(2nh) 1/2 "\ h 

where p G (R) satisfies \\p\\ L 2^ = 1. Then the functions Uh are eigenfunctions of 
H(hD x ) and the measures \uh(xi, X2)\ 2 dxi dx-i obviously concentrate on the diagonal {x\ = 
x 2 }. 

Note that statement (2) of Theorem 11.21 has already been proved in the case H(£) = |^| 2 
and Th = 1/h in [6] and [3] with different proofs. However, the extension of the methods in 
these references to more general H is not straightforward, even in the case where H(£) = 
£-A£, where A is a symmetric linear map : (R d )* — > R d (i.e. the Hessian of H is constant), 
the difficulty arising when A has irrational coefficients. 
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The proof in |3] extends to the (t, x)-dependent Hamiltonian |£| 2 + /i 2 V (t, x) with V con- 
tinuous except for points forming a set of zero-measure. Recently, this has been extended in 
|7J to more general perturbations of the Laplacian (allowing for potentials V G L°°(R x T d )) 
by means of an abstract argument that uses the result in [6j [3] for V = as a black-box. 
In fact, the proof of the result in [7] applies to our context. 

Remark 1.5. Theorem \1.2\ and [7J imply that statement (2) of Theorem \1.2\ also holds for 
sequences of solutions to the Schrddinger equation corresponding to the perturbed Hamilton- 
ian H(hD x ) + hr^ l V{t) where V G L°°(R; C{L 2 {T d ))). The size hr^ 1 of the perturbation 
is in some sense optimal; in Section^ 3 we present an example communicated to us by 
J. Wunsch showing that absolute continuity of the elements of M. {1/h) may fail in the 
presence of a subprincipal symbol of order hP with (3 G (0, 2) even in the case H (£) = |£| 2 . 

Theorem ll.2( 2) admits a microlocal refinement, which allows us to deal with more general 
Hamiltonians H whose Hessian is not necessarily definite at every £ G M. d . Given fi G A4 (r) 
we shall denote by ft the image of fi under the map tt 2 : (x,£) 1 — > £■ It is shown in the 
appendix that fx does not depend on t (it can be obtained as jx = (^2)^1^0, where the 
measure liq is an accumulation point in V (T*T d ) of the sequence (w^ ,))■ 

Theorem 1.6. Let li G AA(l/h) and denote by fx^(t,-) the disintegration of fi(t,-) with 
respect to the variable!;, i.e. for every 9 G L l (R) and every bounded measurable function f: 

[ 0(t) [ f(x,Z)fi(t,dx,dt)dt= fe(t)f ( [ f(x,Z)fit(t,dx))ji(dZ)dL 

Then for fi-almost every £ where d 2 H(^) is definite, the measure fi^(t, •) is absolutely con- 
tinuous. 

Let us introduce the closed set 

Ch '■— {£ G R d : d 2 H(i) is not definite} . 

The following consequence of Theorem 11.61 provides a refinement on Theorem 11.2( 2). in 
which the global hypothesis on the Hessian of H is replaced by a hypothesis on the sequence 
of initial data. 

Corollary 1.7. Suppose v G A4(l/h) is obtained through an h-oscillating sequence (uh) 
having a semiclassical measure liq such that liq (T d x Ch) = 0. Then v is absolutely 
continuous with respect to dtdx. 

We show in Section Hl2 that absolute continuity may fail for the elements of A4 {1/h) 
when H (£) = |£| , k G N and k > 1; a situation where the Hessian is degenerate at £ = 0. 

1.4. Second-microlocal structure of the semiclassical measures. Theorem 11.61 is a 
consequence of a more detailed result on the structure of the elements of A4 {1/h). We 
follow here the strategy of [3J that we adapt to a general Hamiltonian H(£). The proof 
relies on a decomposition of the measure associated with the primitive submodules of (Z d )*. 
Before stating it, we must introduce some notation. 
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Recall that (R d )* is the duafl of R d . We will denote by (Z d )* the lattice in (R d )* defined 
by (Z d )* = {£ G (R d )*,£.n G Z, Vn G Z d }. We call a submodule A C (Z d )* primitive if 
(A) fl (Z d )* = A (here (A) denotes the linear subspace of (R d )* spanned by A). Given such 
a submodule we define: 

(7) J A := {£ G (R d f :dH(£)-k = 0, VA; G A} . 

We note that I\ \ Ch is a smooth submanifold. 

We define also L p (T d , A) for p G [1, oo] to be the subspace of L p (T d ) consisting of the 
functions u such that u(k) = if k G (Z d )* \ A (u(k) stand for the Fourier coefficients 
of u). Given a G (T*T d ) and £ G denote by (a) A (■,£) the orthogonal projection of 
a(-,0onL 2 (T d ,A): 

(8) < a >A = £^(07|^ 

We denote by ^{ a > A (£) the operator acting on L 2 (T d , A) by multiplication by (a) A (•,£). 

Theorem 1.8. Let p G .M (l//i). For every primitive submodule A C (Z d )* i/iere exzsi 
a positive measure p^ £ C (WL;Ai+ (T*T d )) supported on T d x 1^ and invariant by the 
Hamiltonian flow (p s such that : for every a G C c °° (T*T d ) that vanish es on T d x Ch and 
every 6 G L 1 (R) : 

(9) 0(t) a(x,£)p(t,dx,d£)dt = V (t) a(x,g) p A (t,dx,d£) dt, 

JR J T*T d ACZ d M JT d xI A 

the sum being taken over all primitive submodules of (Z d )*. 

In addition, there is a measure p\ on I\, taking values on the set of non-negative, 
symmetric, trace-class operators acting on L 2 (T d , A), such that the following holds: 
(10) 

/ a (x, ha (t, dx, d£) = [ Tr (e~f ^^'^m^ (a) e^ H ^ D ^p K {da) 

JT d xI A Jl A ^ 

When the Hessian of H is definite, formulae (TJ|), / TTOj) /ioZd /or ei>en/ a G C c °° (T*T d ) and 
therefore completely describe p. 

Theorem 11.81 has been proved for H (£) = |£| 2 in [2 1 J for d = 2 and in |3j for arbitrary 
dimension (there, a formula similar to (TTO]) is proved for the x-dependent Hamiltonian 

iei 2 + ^(x)). 

We see that Theorem 11.81 allows to describe the dependence of p on the parameter t. 
As was noticed in |20l [2Tj . the semiclassical measures of the sequence of initial data (uh) 
do not determine uniquely the time dependent semiclassical measure p. On the other 
hand, these are fully determined by the "measures" p\, which are two-microlocal objects 
determined by the initial data (uh)- The contents of f fTUj) is that the measure p\ {t, dx, d£) 
can be obtained by transporting p\ by a certain Schrodinger flow, and then tracing out 



4 Later in the paper, we will tend to identify both by working in the canonical bases of . 
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certain directions. The p\ are obtained by a process of successive two-microlocalizations 
along nested sequences of submanifolds in frequency space; this process gives an explicit 
construction of p in terms of the initial data. This two-microlocal construction is in the 
spirit of that done in [291 fT2| IT3] in Euclidean space. We also refer the reader to the 
articles [33j EH [35] for related work regarding the study of the wave-front set of solutions 
to semiclassical integrable systems. 

Remark 1.9. The arguments in Section 6.1 of [3] show that Theorem \l.b\ is a consequence 
of Theorem II. 8i Therefore, in this article only the proof of Theorem \1.S\ will be presented. 

Remark 1.10. Theorem M. 8\ holds for the time scale = 1/h. 7/r/jC the elements 
of A4 (r) can be described by a similar result (see Section \3. 3\) involving expression (TJ|). 
However, the propagation law appearing in the formula replacing / TlOj) involves classical 
transport rather than propagation along a Schrddinger flow, and as a result Theorem M .2( 2) 
does not hold for Th <C 1/h. 

When the Hessian of H is constant Theorem 11.81 gives a complement to the results 
announced in [3] (where the argument was only valid when the Hessian has rational coef- 
ficients). The statement is as follows. 

Corollary 1.11. Suppose H (£) = £• At; where A : (M d )* — > M. d is a symmetric definite 
linear map. Given v G Ai {1/h) there exists, for each primitive module A C (Z d )*, a non- 
negative, self-adjoint, trace-class operator £a acting on L 2 (T d , A) such that, forb G C (T d ) 
and 6 G L 1 (R); 

(11) fe(t) [ b{x)u{t,dx)dt= V f 0(OTr(m (6)A e-^^S A e itH ( o ^)^. 

In fact, Sa := p\ (Ia), where p\ is given by Theorem M . 81 

Comparing the special case (II ip to the general case (Q, (fTUj) . we see that in the former 
case the propagation law involve the constant propagator e ltH( * Dx \ whereas in the latter 
case we need a "superposition" of propagators e~ 1 ^ d2H ^ Dy ' Dv depending on a G Ia- 

1.5. Hierarchy of time scales. In this section, we present a more detailed discussion 
on the dependence of the set Ai (r) on the time scale r; we shall also clarify the link 
between the time-dependent Wigner distributions and those associated with eigenfunctions. 
Eigenf unctions are the most commonly studied objects in the field of quantum chaos, 
however, we shall see that they do not necessarily give full information about the time- 
dependent Wigner distributions. 

A particular case of our problem is when the initial data (uh) are eigenfunctions of 
H{hD x ). We note that the spectrum of H (hD x ) coincides with H f/iZ d J; given Eh G 
H (hZ d ) the corresponding normalised eigenfunctions are of the form: 

(12) u h {x)= c " e ^ with EKf = 7Ai- 

H(hk)=E h kdlA ^ > 
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In addition, one has: 

Vh, (r h t, •) = \S T f^u h \ = \u h \ 2 , 
independently of (77J and t. Let us denote by Al (00) the set of accumulation points 
in V (T^) of sequences \uh\ 2 where (v,h) varies among all possible /Voscillating sequences 
of normalised eigenfunctions ( JT2|) . Denote by Al av (r) the subset of V (T d ) consisting of 
measures of the formJl 

1 

v (t, •) dt, where v G Conv Al (r). 

Proposition 1.12. Suppose (77J and (t£) are time scales tending to infinity and such that 
r' h r h . Then: 

M (00) C M (r) C L°° (R; Ai av (r')) . 

Remark 1.13. v4s a consequence of Theorem \TJ3( 2) we obtain that all eigenfunction limits 
Al (00) are absolutely continuous under the definiteness assumption on the Hessian of H. 

A time scale of special importance is the one related to the minimal spacing of eigenval- 
ues : define 

(13) t» := h sup [\El- El]' 1 : E l h ± E 2 , E\, E 2 h G H (hZ d )} . 

It is possible to have = 00: for instance, if H (£) = |£| a with < a < 1 or H (£) = £-^4£ 
with A a real symmetric matrix that is not proportional to a matrix with rational entriesjj 
in some other situations, such as H (£) = |£| Q with a > 1, (ITS]) is finite : = h l ~ a . 

Proposition 1.14. If Th ^> one has: 

M (r) = Conv M (00). 

This result is a consequence of the more general results presented in Section 
Note that Proposition 11.141 allows to complete the description of Al (r) in the case 
H (0 = |£| as the time scale varies. 

Remark 1.15. Suppose H (£) = |£| 2 ; or more generally, that ~ and i/ie Hessian 
of H is definite. Then: 

!/ta<1//i, 3z/ G Al (t) sitc/i t/iat z/ _L czMx; 

ifr h ~l/h, M{t) CL^RjL 1 ^)); 

ifr h ^>l/h M (r) = Conv Al (00). 

Note that in this case the regularity of semiclassical measures can be precised. The 
elements in Al (00) are trigonometric polynomials when d = 2, as shown in |18| ; and in 
general they are more regular than merely absolutely continuous, see [U [181 [28]. The same 
phenomenon occurs with those elements in Al (1/h) that are obtained through sequences 
whose corresponding semiclassical measures do not charge {£ = 0}, see [2]. 

5 ConvAT stands for the closed convex hull of a set X C L°° (WL;V (T rf )) with respect to the weak-* 
topology. 

6 This is the content of the Oppenheim conjecture, settled by Margulis [51 123|. 
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1.6. Organisation of the paper. The key argument of this article is a second microlo- 
calisation on primitive submodules which is the subject of Section 12.21 and leads to Theo- 
rems 12.41 and 12.51 In Section |3l successive microlocalisations allow to prove Theorem 11.81 
and Theorem 11.2( 2) when ~ 1/h. Examples are developed in Section H] in order to 
prove Theorem 11.2( 1). Finally, the results concerning hierarchy of time-scales are proved 
in Section [5] (and lead to Theorem 11.21 for ^> 1/h). 

Acknowledgements. The authors would like to thank Jared Wunsch for communicating 
to them the construction in example (3) in Section Hl3. They are also grateful to Luc 
Hillairet for helpful discussions related to some of the results in Section Part of this 
work was done as F. Macia was visiting the Laboratoire d'Analyse et de Mathematiques 
Appliquees at Universite de Paris-Est, Creteil. He wishes to thank this institution for its 
support and hospitality. 

2. TWO-MICROLOCAL ANALYSIS OF INTEGRABLE SYSTEMS ON T d 

2.1. Invariant measures and a resonant partition of phase-space. As in [3J, the 
first step in our strategy to characterise the elements in Ai (r) consists in introducing a 
partition of phase-space T*T d according to the order of "resonance" of £, that induces a 
decomposition of the measures fi £ A4 (r). We say that a measure fi £ M+ (T*T d ) is 
a positive H -invariant measure on T*T d whenever \i is invariant under the action of the 
hamiltonian flow 

(14) (0 s )^ = /i, with cj) s (x, = (x + sdff(0,0 ■ 

Recall that £ is the family of all primitive submodules of (Z d )* and that with each A £ C, 
we associate the set I\ defined in (J7|): if A -1 C M> d is the orthogonal to A with respect to 
the duality in (R d )* x R d then J A = dH~ l (A L ). Denote by Qj C R d , for j = 0, ...,d, the 
set of resonant vectors of order exactly j, that is: 

n d := {££(ffi d )*:rkA 5 = d-j}, 

where 

A 5 := {k £ (Z d )* : k-alH^) = 0} . 

Note that the sets flj form a partition of (IR d )*, and that Q = dH^ 1 ({0}); more generally, 
£ £ Qj if and only if the Hamiltonian orbit {<p s (x, £) : s £ M} issued from any x £ T d in 
the direction £ is dense in a subtorus of T d of dimension j. 

The set Q := Uj=o ^? ^ s usuan y called the set of resonant momenta, whereas fl^ = 
(M. d )* \ Q is referred to as the set of non-resonant momenta. 

Finally, write 

Saying that £ £ is equivalent to any of the following statements: 

(i) for any x £ T d the time-average |; J Q T S Xo+t dH(£) (%) dt converges weakly, as T — > oo, 
to the Haar measure on the torus Xo + T A ±. Here, we have used the notation 
T A x := A L / (2nZ d n A 1 ); 
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(ii) A 5 = A. 

Moreover, if rkA = d - 1 then R A = dH' 1 (A- 1 \ {0}) = I A \ Q . Note that, 
(15) (R d T = U «A, 

that is, the sets i?A form a partition of (R d )*. As a consequence, any measure /i G 
.M-^T*!^) decomposes as 



(!6) /^X^t'x-Ra- 



Therefore, the analysis of a measure /i, reduces to that of ^\t<i xRa for all primitive submod- 
ule A. Given an H- invariant measure /i, it turns out that //|T d xR A are utterly determined 
by the Fourier coefficients of /i. Indeed, define the complex measures on M. d : 

p ^—ik-x 

fik-= - — -jj^nidx,-), keZ d , 

JT d (27TJ ' 

so that, in the sense of distributions, 

ik-x 



kei d 

Then, the following Proposition holds. 



Proposition 2.1. Let // G A4 + (T*T d ) and A G £. 77&e distribution: 



Ak-x 



d/2 



is a finite, positive Radon measure on T*T d . 

Moreover, if n is a positive H -invariant measure on T*T d , then every term in the decom- 
position l[Tb)) is a positive H -invariant measure, and 

(17) tAj d xR A — (^aI^xRa- 

Besides, identity j[T7\) is equivalent to the fact that n1fd xRA is invariant by the translations 

(x, £) i — > (x + v, £) , for every v G A x . 

The proof of Proposition 12.11 follows the lines of those of Lemmas 6 and 7 of [3]. We 
also point out that this decomposition depends on the function H through the definition 
of I\. Consequently, it is possible to make different choices for decomposing a measure /i; 
this fact will be exploited in Section [3j 
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2.2. Second microlocalization on a resonant submanifold. Let (uh) be a bounded 
sequence in L 2 (T d ) and suppose (after extraction of a subsequence) that its Wigner dis- 
tributions Whif) = w h tTh converge to a semiclassical measure /x G L°° (R; M. + (T*T d )) in 

h u h 

the weak-* topology of L°° (R; V (T*T d )). 

From now on, we shall assume that the time scale (77J satisfies: 

(18) (hTh) is a bounded sequence. 

Given A G £, the purpose of this section is to study the measure H~\r d xR A by performing a 
second microlocalization along I\ in the spirit of [T2| [T3| [HI |29l |27J and [3] [21]. By Propo- 
sition ETTJ it suffices to characterize the action of ^\t<i xRa on test functions having only 
x-Fourier modes in A. With this in mind, we shall introduce two auxiliary "distributions" 
which describe more precisely how Wh (t) concentrates along T d x I\. They are actually 
not mere distributions, but lie in the dual of the class of symbols S\ that we define below. 

In what follows, we fix £ G R\ such that d 2 H(^ ) is definite and, without loss of 
generality, we restrict our discussion to normalised sequences of initial data (uh) that 
satisfy: 

u£(Jfe) = 0, for hk G m. d \B(£ ;e/2), 

where B(^ ,e/2) is the ball of radius e/2 centered at £o- The parameter e > is taken 
small enough, in order that 

d 2 H(i) is definite for all f G B(£ Q , e); 

this implies that /aHB((o, e) is a submanifold of dimension d — rkA, everywhere transverse 
to (A), the vector subspace of (R d )* generated by AH By eventually reducing e, we have 

B(£„,e/2) C (/ A n%£))®(A), 

by which we mean that any element £ G B(^ , e/2) can be decomposed in a unique way as 
£ = a + 77 with a G /a H -B(£ , e ) an d ?7 G (A). We thus get a map 

(19) F : B(£ , e/2) — ► (J A n B(£ , e)) x (A) 

£— »j(0) 

With this decomposition of the space of frequencies, we associate two-microlocal test- 
symbols. We denote by S\ the class of smooth functions a (x, £, rj) on T*T d x (A) that 
are: 

(i) compactly supported on (zc, £) G T*T d , £ G B(£o,e/2), 

(ii) homogeneous of degree zero at infinity w.r.t. 77 G (A), £e. such that there exist 
R > and a hom G C c °° (T*T d x S(A)) with 

a (x, £, 77) = a hom ( x, £, |-r ) , for > i? and (x, £) G T*T d 



7 This can be made by applying a cut-off in frequencies to the data. 

8 Note that this is achieved under the weaker hypothesis that d 2 H(^) is non-singular and defines a 
definite biliear form on (A) x (A) . 
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(we have denoted by S(A) the unit sphere in (A) C 

(iii) such that their non vanishing Fourier coefficients (in the x variable) correspond to 
frequencies k G A: 

keA 



(2ir) d/2 ' 



We will also express this fact by saying that a has only x-Fourier modes in A. 

Let x £ C-T ((A)) be a nonnegative cut-off function that is identically equal to one near 
the origin. For a G S\, R > 1, 5 < 1, we decompose a into: a(x,^,i]) = ^| =1 (a:, ^, 77) 
with 

n Jr £ „\ .— n (r C „\ {l _ ^ \ ( 1 _ ^ ( T ±_ 



ai (x,Z,ri) := a(x,e,77)(l-x(|)) (l - 

(20) a 2 (x,t, V ) := a(x, Z, v ) (l - x (|)) X (^) , 

(21) a 3 (x,^,r]) := a(x, f, r])x ■ 

This induces a decomposition of the Wigner distribution: 

Mt) = w T h \ s (t) + w w (t) + w'fag (t) 
(when testing against functions a with Fourier modes in A), where: 

{ W W (t) > a ) '■= / «2 (x, f , T h r)(£)) w h (t) (dx, d£) , 

JT*T d 

(22) (wi A , h ,R (t) , a) := / a 3 (x, f, r^(0) w h (t) (dx, df) , 

JT*7 d 

and 

(23) \ w h\s (*) > a ) := / fl i ( x > WO) w fc (*) ( dx > d > 

that we shall analyse in the limits ft, — >■ + , i? — >■ +00 and 5 — > (taken in that order). 
One sees that 



lim lim lim f 9(t) (w[ K RS (t) ,a) dt = [ [ 0(i)aoo ( x, £, t^tt ) fi(t, dx, d£)] T d XI cdt 



where |i6M (77J is the semiclassical measure obtained through the sequence (uh). The 
restriction of the measure thus obtained to T d x R\ vanishes, and we do not need to further 
analyse the term involving the distribution w h A R5 (t). 
For a G S\, we introduce the notation 

Op£(a(a:, f , 77)) := Op h (a (x, f , r^(0)) 
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so that the distributions wjf RS (t) and Wi Ai h,n{t) can be expressed for all t G M by 



iw 



(t),a) = {u^S^O^a^Sl^Uh)^^). 



Notice that, for all f3 G N d , 

d^(a(x,h^r hV (hO)) <C p {r h hf\. 

S Loo 

The Calderon-Vaillancourt theorem (see |8j or the appendix of [3] for a precise statement) 
therefore ensures that there exist N G N and Cjy > such that 

(24) Va G Si, ||0p£(a)|| £(i2(Rd)) < C N ^ ||^a|U», 

|a|<JV 

since (h,Th) is bounded. 

As a consequence of (I24p . both u^ A r<5 anc ^ w IaAR are bounded in L°° (R; (<S A ) ). After 
possibly extracting subsequences, we have for every tp G L 1 (R) and a G S\, 

ip (t) </i A (t, •) , a) dt := lim lim lim / <p (t) (w^ R 5 (t) , a) dt, 

and 

(25) / ip(t)(fj, A (t,-),a)dt:= lim lim / ip (t) (w lAAR (t) ,a) dt. 

Remark 2.2. When <C /j/ie quantization of our symbols generates a semi-classical 
pseudodifferential calculus with gain hr^. The operators Op^(a) are semiclassical both in 
£ and 77. 27ms implies that the accumulation points ft A and are positive measures (see 
for instance |14] ). 

Because of the existence of Rq > and of ahom G (T*T d x § (A)) such that 

a (x, f , r/) = a h0 m £, J , for |r/| > R , 

for large enough, the value (w^ RS (t) , a) only depends on ahom- Therefore, the limiting 
object fi K {t, ■) G (Si)' is zero- homogeneous in the last variable 77 G M d , supported at 
infinity, and, by construction, it is supported on £ G /a- This can be also expressed as the 
fact that jl k is a "distribution" on T d x J A x (A) (where (A) is the compactification of (A) 
by adding the sphere §(A) at infinity) supported on {77 G 8(A)}. Besides, the distribution 
fl\ is supported on T d x x (A). Indeed, we have for all t, 

(26) (w lAth>R (t), a(x, f , 77)) = (w lAAR (t), a(x, <r(f ), »?)) + ^(r^ 1 ) 
since, by (124|) . 

Op£(a 3 (x,£,77)) = Opt(a(x,a(0+T^,77)x(W^)) 

= Op^(a(x,a(0,r/)x(V J R)) + O(r, 1 ) 
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where the O(t)^ 1 ) term is understood in the sense of the operator norm of £(L 2 (M d )) and 
depends on R (the fact that we first let h go to + is crucial here). 

From the decomposition Wh{t) = w h A R s (t) + Wi At h,n (t) + w h A R s (t) (when testing against 
symbols having Fourier modes in A), it is immediate that the measure /i (t, -)lT d xi?A * s 
related to fi A and jl^ according to the following Proposition. 

Proposition 2.3. Let 

fi A (t, ) ■= I fl A {t,-,dr])] T d xRA , /i A (t, •):=/ ft A {t,-,dr])] T d xRA . 

J (A) J (A) 

Then both /i A (t, •) and /iA (t, ■) are H -invariant positive measures on T*T d and satisfy: 

(27) A*(V)1t*x* a = A* A (*,■)+ A*a (*,■)■ 

This proposition motivates the analysis of the structure of the accumulation points 
/iA (t, ■) and ft A (t, •). It turns out that both fi A and /i A have some extra regularity in the 
variable x, although for two different reasons. Our next two results form one of the key 
steps towards the proof of Theorem 11.21 

Let us first deal with ft A (t, •). We define, for (x,£,rj) G T*T d x ((A) \ {0}) and s£l, 

0° {x,C,v) ■= (x + sdH(£U,ri), 



<t>\ (x, e, V) ■= U + sd 2 H(a(0) ■ |^| , C, f]j 



This second definition extends in an obvious way to 7/ G S(A) (the sphere at infinity). On 
the other hand, the map (x, £, if) i— > <p\, (x, £, rj) extends to r\ — 0. 

Theorem 2.4. jx A (t, •) is a positive measure on T x I\ x (A) supported on the sphere at 
infinity §(A) m £/ie variable rj. Besides, for a.e. t G R, £/te measure fi A (t,-) satisfies the 
invariance properties: 

(28) (0°)> A (*,•) = //(*,•), (#),/i A (t,-) = /i A (V), sei 

Note that this result holds whenever <C l//i or r/, = l//t. This is in contrast with 
the situation we encounter when dealing with jj\(t, •). The regularity of this object indeed 
depends on the properties of the scale. 

Theorem 2.5. (1) The distributions fi^{t,-) are supported on T d x I\ x (A) and are 
continuous with respect to t G M. Moreover, they satisfy the following propagation law: 

(29) Vt G R, £ A (t, £, r/) = (0^)^(0, x, £, 7?). 

(%) // Th <C l//i t/ien /tA (t, •) is a positive measure. When Th = 1/h, the projection of 
/iA (t, ■) on T*T d is a positive measure, whose projection on T d is absolutely continuous 
with respect to the Lebesgue measure. 

Remark 2.6. For Th = 1/h the propagation law satisfied by distributions /i A (t, •) can be 
interpreted in terms of a Schrodinger flow type propagator. The precise statement can be 
found in Proposition \2.13\ in SectionlKJ. 
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Remark 2.7. Note that for all £ £ (R d )* \ Cjj (recall that Ch stands for the points where 
the Hessian d 2 H (£) is not definite) we have R d = A 1 - © d 2 H (£) (A). Therefore, the flows 
4>° s and <fil are independent on T d x (R\ \ Ch) X (A). 

Remark 2.8. If rk A = 1 then < f^5j) implies that, for a.e. t £ R ; and /or any z/ £ (A), tae 
measure jl A (t, •)lT d xi? A x{A) ^ s invariant under 

(x, a, rj) i — > (x + d 2 H(a) ■ is, a, rj). 

On the other hand, the invariance by the Hamiltonian flow and Proposition ^. 1\ imply that 
jl A (t, -)lT d x_R A x(A) is a l so invariant under 

(x,a,rj) i — > (x + v,a,rj) 

for every v £ A -1 . Using Remark \2. 71 and the fact that the Hessian d 2 H(o~) is definite 
on the support of fi h (t, •)lT d x_R A x(A); we conclude that the measure jl K (t, -)]T d xR A x{A) 
constant in x £ T d in this case. 

Remark 2.9. Consider the decomposition \i (t, •) = YIa&c^ ') + ^Ae£/ iA ') ■ 9^ ven 
by Proposition \2.3l When = 1/h, Theorem \2.5\ implies that the second term defines 
a positive measure whose projection on T d is absolutely continuous with respect to the 
Lebesgue measure. 

We now give the proof of Theorems 12.41 and 12.51 

2.3. Invariance properties of jl A . In this section, we prove Theorem 12.41 The positivity 
of fi k (t, •) can be deduced following the lines of [14J §2.1, or those of the proof of Theorem 1 
in [16]; see also the appendix of [3]. 

Let us now check the invariance property (1281) . We use the following Lemma which gives 
approximate transport equations by the flow (fP s . 

Lemma 2.10. For every a £ S\ and (p £ Cg°(R) ; we have 

/ ip(t){u h ,S^ Tht Op£(a)5£** u h ) L 2 {Jd) dt = / ip{t)(u h , Op^(ao^° Tht )u h ) L 2 {Jd) dt + o(l). 
Jr Jr 

Remark 2.11. Consider the h- dependent flow 

<f>l' h {x, f , n) := (x + sr h (dH(£) - dH(a(£))) , £, rj) , 
then if a has only Fourier modes in A ; we have 

a o (j) 1 / (x, £, rj) = MZ> r})e ik < x+STh ( dH ®- dH ^))) 
fceA 

= v)e ik - (x+SThdm)) = a o ^ rhS (x, £, rj) . 

fceA 

This comes from the fact that for every A; £ A and £ £ M. d , one has k ■ dH = 0. 
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We postpone the proof of Lemma r2.10l to the end of this section and we start proving (|28|) . 
The invariance by the "geodesic flow" 0° is standard and can be proved following the lines 
of the proof of property 3 in the appendix. Using (120|) . we have 

(30) / <p(t)(jl A (t, •), a)dt = lim lim lim f tp(t) (w'fa (t) , a) dt 



= lim lim lim / (p(t)(w h (t),a 2 (x,£„T h rj (£,))) dt 

(along subsequences). Notice that the symbol 

02 o 4>l{x, £, rj) = a 2 (^x + sd 2 H(a(i))^ £, r^j , 

is a well-defined element of S\, since, for fixed R, a 2 is identically equal to zero near r] = 0; 
moreover 

Vw G §(A), (a 2 o 4>l) hom (x, £, u) = a hom (x + sd 2 H(a(£))u, £, u); 

therefore, 

(31) / V (t)(ji A (t, •), a o (j>\)dt = lim lim lim / <p(t) (wi\ 5 (f),oo 0, 1 ) dt. 

In order to relate fl3T|) to f l30|) we note that the symbol: 

a 2 o <t> l * H (x, £, 77) , 

satisfies: 

rlj< 

I r > ; n = i7 .-, I / ' ■ ■ 

WO I 



«2 o (t) L f/\ ri \ (x, f, rj) = a 2 [x + y-T—(dH(£)-dH(a(£)),£,T h ri(£,) 



(32) = (a 2 o0i)(x,e,r^(O) + O((5); 

we have used that, on suppa 2 , we have < C5 and the function 

(33) G(0 = / d 3 H(a(£) + tr){0)(l - t)dt, 

Jo 

is uniformly bounded. On the other hand, by Lemma 12.101 and Remark I2.11[ we have 

(34) (wffag (t),a) = (w h (t),a 2 (x,Z,T h r)(Z))) = (w h (0) , a 2 o <j>\> h (x, £, r^))) + o h (l). 

Therefore, combining (132]) and (I34p we obtain: 

(35) 

(vfau (t),ao (Pi) dt= [ <p(t) (w h (0) , a 2 o 0i o $ h (x, £, r hV (£)))dt + 0(5) + o h {l) 



: [t-^] R(0),a 2 o <ft h {o:i,. n ,/(4)))r/f ■■ ■()!<)!■ ■■ <>,,(!). 
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Since \r]\ > R on the support of a 2 , we have for all K e N 



lim 

R— >+oo 



sup sup 

x,€,V at€N 3d ,\a\<K 



d 



{<p{t - s/\r)\) - tp(t)) (a 2 o t 1,fc ) , 77) 



<it = 0, 



which implies, in view of fl35l) : 

lim I™ lim / (*) , a ^> ^ 

<5— >0 R—>-+oo h—tO Jjj 

= lim lim lim / (p(t)(w h (0) , a 2 o <f)]' h (x, £, T h r}(£)))dt. 

Applying again ( l34|) to the left hand side of the above identity concludes the proof of 
Theorem I2.4I 

Let us now prove Lemma I2.10I 

Proof of Lemma \2.KA Write 
(36) 

<p{t){u h ,S- Tht Op£(a)S^ u h ) L2{Td) dt 

H(hk) -H(hj) 



k,j£Z d ,k-jeA 



h 



Uh (k) u h (j) aj-k h——, r h r) h 



k+j 



k + j 



where a(x, £,77) = X^eA^ fc (£' r ?) e ■ Notice that, 



H(hk)-H(hj) = hdH [h 



k + j 



(k-j) + h d r h [h 



k + j 



k-j 



where satisfies, for every .K" C M. d compact and convex and for every /3 G N d 



3f r fc (£, 



< C 



K,/3 



f , hi e AT. 



Therefore, since (p is uniformly Lipschitz, we have, for hk, hj G K 



H (hk) - H (hj) 



h 



(p r h dH h 



k + j 



(k-j)\+h 2 T h M h [h—^ } k~j 



k + J 



with 

M h (Z,£) =r h (Z,£) [ $ (r h dH(0-i + sh 2 r h r h (^£))ds. 
Jo 

Plugging this equation in ( |36|) we obtain: 

ip(t)(u h ,S^ Tht Op£(a)S^* u h } L 2 {Td) dt 

<p(t) (u h , Op£ (a (x + T h tdH(0,£, r})) dt + h 2 r h (u h , Op£(R h a ) 
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where the symbol R%(x,£,r)) is characterized by its Fourier coefficients : 

Since the function a is compactly supported in (x, £), we deduce from the properties of 
and rh that for all (3 £ N d , there exists Cp > such that 



Passing to the limit h — > + allows to conclude the proof of the Lemma. □ 

2.4. Propagation and regularity of fi^. This section is devoted to proving Theorem 12.51 
and showing that, when = 1/h, the distribution fl\ satisfies a propagation law that 
involves the family of unitary propagators 

We start by proving Theorem 12.5( 1). The statement on the support of jl^ was already 
proved in Section 12.21 The propagation law (and hence, the continuity with respect to t) 
comes from the following result. 

Proposition 2.12. For every a £ S\ and every t £ R the following holds: 

(wi A ,h,R (*) > a ) = \ w IaAR (°) > a ° + °^ C 1 ) > 

where 

(37) 0, 1 (z, f , v) ■= (x + td 2 H{a(C))v, C, v) ■ 
Proof. We use Lemma 12.101 and Remark 12.111 to conclude : 

(wi A ,h,R (*) ' a ) = (°) > a 3 ° <f>t' h ( x ,£,i~ h r)(£))^ + o ft (l), 
where a 3 is defined by (|2ip . By definition of </>J ,?i and by Taylor expansion, we obtain 
a 3 o (x, f , T h ri{0) = a 3 (x + tr h (dH(C) - dH r hV (0) 

= a 3 (x + td 2 H(a(0)r hV (0 + tr h G{£) [n(S)MS)U, ^(0) 

=: 6 ft (t,X,£,T h T](£,)) , 

where G is defined by (1531 and is bounded and smooth on the support of a 3 . Therefore, 

(38) 6ft(t, x, £, tj) = a 3 (x + td 2 H(a(0)v + r h l tG(0 [ V , rj\,£, rj) 

= a 3 (x + td 2 H(a(0H £, 7?) + O^- 1 ), 
from which the result follows. □ 

Let us now focus on statement (2) of Theorem 12.51 The result concerning time scales 
Th ^ 1/h was already discussed in Remark 12.21 

From now on, suppose = 1/h. Let us introduce some notations. We consider the 
set L 2 (T d ,A) C L 2 (T d ) consisting of functions having only Fourier modes in A. For any 
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a G S A , and for any a G I\, we can associate the operator Op 1 (a a (y,r])) which acts on 
L 2 (T d ,A), defined as follows. For $ G L 2 (T d ,A), 

Vl (a a (y, V ))$ = ^ + $ (") e^ +A)i/ , 

where a(x,£,T]) = XlfcgA ®k (£' V) (J^ws ■ Thus 0p 1 (a o .(y, 77)) is the Weyl quantization of the 
symbol a a (y,rj) := a(y,a,rj). In particular, for a = a(x, £) independent of the variable 77, 
OPi ( a o-(y)) is the multiplication operator 

(39) Op 1 My))$ = a(y,<r)$ 

(we will simply denote by a a (y) this multiplication operator). We have the expression 
similar to formula (|52p in the appendix, 

^ ^ / y _|_ y/ \ ^ — 

($,0 Pl (a ff (y,r7))$) 

K— ) $ fa) • 

^ ^ u,u'eA ^ ' 

Then, the last statement of Theorem 12.5( 2) is a consequence of the following Proposition^ 

Proposition 2.13. There exists M G M+ (l A ; C 1 (L 2 (T d , A))) stic/j that for all a G 
C£°(T*T d ) iwt/i Fourier modes in A and all if G L l (M), 

[ if (t) (/2 A (t, •) , fl)dt = / ^(t) / Tr ( e -^H{a)D v Dy^y^H^DyD^^A ^ 

Jm Jr Ji a ^ ' 

Remark 2.14. (i) The operator-valued measure M is globally defined, it describes the limit 
of (wi A) h,ii(t), a) for symbols a = a(x,£). For symbols a = a(x, £, 77) G S\, one cannot build 
such a global measure (see Remark \2.19i which emphasizes the technical obstruction), 
(ii) Note that for any given a, the operator e^ d2H ^ Dy Dy obviously preserves L 2 (T d ,A). 

The proof of Proposition 12.131 relies on three steps: 

(1) We first define an operator Uh which maps (27rZ d )-periodic functions on (27rZ d )- 
periodic functions with Fourier frequencies only in A. 

(2) Then, we express Wj A: h,ii(t) in terms of Uh and the operators e^ d2H ^ Dy ' Dy . 

(3) We then conclude by passing to the limit when h —> and R — > +00. 

First Step: Construction of the operator Uh- We introduce an auxiliary lattice A C Z d 
such that A ± © A = Z d (recall that A^ is the orthogonal of A in the duality sense). We 
denote by a the projection on (A), in the direction of A -1 . We have a(Z d ) = A C Z rf . 
For a G (R d )*, we shall denote by cr Q G (A) the linear form cr a (y) = o ■ a(y). We fix a 
bounded fundamental domain D\ for the action of A on (A). For 77 G (A), there is a unique 
M G D A (the "fractional part" of 77) such that 77 - {77} G A. Finally, take b G C^((R d )*) 

9 Recall that given a Hilbert space H, C 1 (H) stands for the space of bounded trace-class operators 
acting on H and, for a Polish space X, M. + (A"; C 1 (H)) denotes the set of positive measures taking values 
on C 1 (H). 
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supported in the ball B(£ , e) C (R d )*, and identically equal to 1 on B(£ ,e/2). We set for 
/ G L 2 (T d ), a G I A , y G T d , 

U h f{a,y) = (2ir)-%e^- y J f(x) + v)^ v ' y e~^ a+v>x dx 

V e(A),(a,r))eF(hZd) 



■ne(A),(cr,T l )€F(hz d ) 



Recall that F is the local coordinate system defined in (Tl9l) . with the property that if 
F (0 = (a (0 , V (0) ^en £ = a (0 + V (0- Note that y) = if (a, 77) £ F(/iZ rf ) 

for every rj G (A) (since the sum has an empty index set). The role of the term e^hr^y 
becomes clear in the following lemma. 

Lemma 2.15. Iff is {2n r £) d -periodic, then Uhf is (27rZ) d -periodic and has only frequencies 
in A. Therefore, Uh maps L 2 (T d ) into the subspace L 2 (T d , A) of L 2 (T d ). 

Proof. It is enough to show that for any a G I\, r\ G (A) such that a + 1] G KL d , 

By definition, l + G (A), and we want to prove that for any k G 27rZ d , (f + {x}) - ^ G 
2irZ. We write 

Hx})K^+{?B)* 

and we know that ■ k G 27rZ. We then use the fact that there exists A G A such that 
{£} = x + A, and write 

ff a ] , cr • k a a , cr , , , . , . . , 

— > ■ A; — = - • a(fc) H-A-A; — — - A;= — • (a(k) -k) + \-k. 

h J h h h h 

Since a+77 = /?i for some / G Z d and since k G (27rZ) a! , we obtain f-(a(/c)— /c) = l-(a(k) — k) 
with a(k) -keA^f) (2vrZ d ) (since a(Z d ) = A C Z d ). Finally, we get f • (a(k) -k)e 2vrZ 
and we also have A • k G 27rZ, which concludes the proof. □ 

Note that if / is (27rZ) d -periodic, the Fourier coefficients of Uhf satisfy 
V V G A, fhf(a, V ) =b(a + hr,-h^Yjf(^ + r,-^j- 

and we have the following Plancherel-type formula. 
Lemma 2.16. If f is (27rZ) Q '-periodic, then we have 

V^G/a, J2\f(k)b(hk)\ 2 = Yl I \Ukf(a,y)\ 2 dy. 

- - •' ae<r(h1 d ) J 1 
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Proof. We have for all (27rZ d )-periodic function /, 

2 



keZ d a€l A ,rie(A) a+r]€hZ, d 



f 



a + rj 
h 



X exp ( - (y (rj - j]')) ) dy 



I \Uhf(c-,y)\ 2 dy. 

ada{KL d ) 1 



□ 

Second step: Link between wi K) h,R and Uh- It is in this step that we really see the relevance 
of the objects introduced previously. It comes from the two following lemmas : 



Lemma 2.17. For any a e S\, 

fl 3 ( x, f , ) w h Uh (dx, d£) 



T*T d 



h J Uh 

{U h u h (a,y), 0p x (a CT (y,v~ {^11 X ( - — ] ] U h u h {a,y)) Lm 



a€a(hZ d ) 

Proof. We have by (T26l) . 

v(0 



R 



+0(h) 



«3 



Then, using (1521) 

03 L a ^rt£L\ w h Uh (dx,do 



w h Uh (dx, d0 = J a 3 U ^ J <(cfc, dfl + O(fc) 



(27T) 



We write 



~~y[/2 u h (k)u h (k')a k ^ k (a (h 



k + k'\ 1 / k + k'\\ ( 1 fk + k' 

— hh v [ h — )) x [hR ri [ h — 



hk = a + hq, hk' = a + hrf , (cr, /177) , (a, /w/) E F(ft,Z a 
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using the fact that k! — k G A. In particular, a (h^^-) = a. Then, 



T*T d 



a 3 (x,a(0,^f)w h Uh (dx,dO 



) , u h ^-+rjju h \j- + r]'Ja rl >- rj la 

rrJ-hW Chid > 



v ' a+hri,(7+ht]'£hZ d 



(2vr) rf / 2 

V 7 ff 6«r(ftZ'') T+ .{^L},^ + {^L}eA 



v + v'\ fv + v' 



which is the desired expression. □ 
To simplify the notation in the computations that follow, we set: 

A(a, V ) := ^d 2 H{a) V - V . 
Lemma 2.18. For any a G S\, for any t G M., 

«3 (x, f , ^P-) w h (t, dx, d£) 



T*T d 



h 



= E (e- iU(CTA) ^(^,y), 

Proof. We use Proposition 12.121 and apply Lemma 12.171 to the symbol bo(t,x, £,,■)]) which 
was defined in the proof of Proposition 12.121 (see fl38|) with h — 0) . Then, the result follows 
by using the fundamental property of the Weyl quantization, namely that, since A (a, if) 
is quadratic in r\ and does not depend on y: 

□ 



Pl (b (t, y, a, ry)) = e~ itA ^ Pl ( a a ( y, V - ^ — M X [ l R h > ) ) e 



24 N. ANANTHARAMAN, C. FERMANIAN, AND F. MACIA 

Third step: Passing to the limit. If a G S\ is compactly supported in 77, the map a 1 — > 
0p 1 (a o .(y, 77)) belongs to the Banach space C C (I\, K(L 2 (T d , A)))0 The dual of this space 
is Ma := M(/a; >C 1 (L 2 (T d , A))), the space of trace-class operator valued measures. Let us 
consider the element p h G .Ma, defined by letting 

(p h ,K)= (U h u h (a,y),K(a)U h u h (a,y)) 

a£cr(hZ d ) 

for all K eC c (l A ; K (L 2 (T d , A))). Lemma 02 implies that (p h ) is bounded in M A if K) 
is bounded in L 2 (T d ). Besides, each p h is positive (meaning that (p h ,K) > if K (a) > 
for all a). We consider M(a) a weak-* limit of the family p h in Ma and we now restrict 
our attention to symbols a G S\ that do not depend on 77. We write 



i? y \RJ R Jo \R R I • & J / I ^ 
and we obtain 

Pl (6o(t, y, a, 77)) = e-^C'AO Pl (a (7 (y) X (r / /i?)) e ^(^») + 0(l/fl). 

whence 



lim lim / a 3 (x, £, — — )wh(t,dx,d£) 

_R->+oo /i->0 Jj>*jd h 



= lim (M(o-), e- i4A(<T '^ Pl KM^/i?)) e**^) . 

Remark 2.19. The arguments used to get rid of the term {^rl crucially exploit the pres- 
ence of a factor l/R in front of it. Such an argument cannot be used for symbols a G 
which depends non trivially of the variable 77 as in Lemma \2.ll\ In such a situation, by 
working in L 2 (M. d ), one can define locally an operator-valued measure M; however, this 
object cannot be globally defined on the torus. 

3. AN ITERATIVE PROCEDURE FOR COMPUTING p 

3.1. First step of the construction. What was done in the previous section can be 
considered as the first step of an iterative procedure that allows to effectively compute 
p(t, ■) solely in terms of the sequence of initial data (uh). Recall that we assumed in §2.2} 
without loss of generality, that the projection on £ of p (t, •) was supported in a ball 
contained in ~R d \ Cjj- We have decomposed this measure as a sum 

p(t,-) = ^/i A (t,-) + X> A (V), 

Aec Aec 

where A runs over the set of primitive sub modules of Z d , and where 

MV)=/ &A(t,-,dr])]rd xRA , p A {t,.)= _p A (t,-,dr])] Jd><RA . 

J (A) J (A) 



10 



Here JC (H) denotes the space of compact operators acting on a Hilbert space H . 
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From Theorem 12. 5[ the distributions jl\ have the following properties : 

(1) /2a(£, dx, d£, drj) is in C (R; (Sjj)') and all its x-Fourier modes are in A; with respect 
to the variable £, /2a ^ 3 drj) is supported in /a; 

(2) if Th -C then for every t G R, /I a (£, •) is a positive measure and: 

/^a (*,•) = £a(0,-), 

where: 

# : (x, f , ?7) 1— ► (ac + sd 2 H{a{i))r 1 , Z, i); 

(3) if r ft = then j {A) jl A {t,; drj) is in C(R; .M + (T*T d )) and ^ dx{A) ji K {t,;d£,drj) is 

an absolutely continuous measure on T d . In fact, with the notations of Section 12^ 
we have, for every a G Cf (T*T d ) with Fourier modes in A, 

/ a(x,Z)ii A {t,dx,dZ,drj)= [ Ti(a a e-^ d2H ^ D y- D yM(da)e^ d2H{(T)D y- D A 

</T d x/ A x{A) Jl A ^ ' 

where M G M.+ C 1 (L 2 (T d , A))) and a CT is the multiplication operator by 
a(-,a), acting on L 2 (Y d , A). 
On the other hand, the measures /i A satisfy: 

(1) for a G iS A , (/i A (t, dx, d£, d?y), a(x, £, 77)) is obtained as the limit of 

«V (*) . «> = jf^ X (^) (l-X ( Z ^) ) « C TW7(0) «fo (*) {dx, dfl , 

in the weak-* topology of L°°(R, («S^) ), as h — > + , R — > +00 and then 5 — > + 
(possibly along subsequences); 

(2) jx A (t, dx, d£, drj) is in L°°(R, M + {T*T d x (A))) and all its x-Fourier modes are in A. 
With respect to the variable 77, the measure ji A (t, dx, d£, drj) is O-homogeneous and 
supported at infinity : we see it as a measure on the sphere at infinity S(A). With 
respect to £ it is supported on {£ G I a}', 

(3) ji A is invariant by the two flows, 

</>° a :(x,Z,7i)^(x + sdH(Z),Z,7i), and # : (x, £, 77) .— >• (x + sd 2 H(a(Z))^-, £, 77). 

This is the first step of an iterative procedure; the next step is to decompose the measure 
/i A (t, ■) according to primitive submodules of A. We need to adapt the discussion of [3]; to 
this aim, we introduce some additional notation. 

Fix a primitive sub module A C Z d and o G I\\ Ch- For A2 C Ai C A primitive 
submodules of (Z d )*, for rj G (Ai), we denote 

A v ((J,Ai) := (Ai®Rd 2 H(a).ri) ± n(Z d Y 
= (Rd 2 H(a).r ] ) ± n A 1; 
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where the orthogonal is always taken in the sense of duality. We note that (a,Ai) is a 
primitive submodule of A 1; and that the inclusion A,, (a, A x ) C A x is strict if rj ^ since 
d 2 H( a) is definite. We define: 

RaM) ■= {^e(Ai),A,( ( x,A 1 ) = A 2 }. 

Because d 2 H(a) is definite, we have the decomposition (R d )* = (d 2 H(a).A 2 ) ± © (A 2 ). We 
define ^° be the projection onto (A 2 ) with respect to this decomposition. 

3.2. Step k of the construction. In the following, we set A = Ai, corresponding to step 
k = 1. We now describe the outcome of our decomposition at step k (k > 1); we will 
indicate in §3.31 how to go from step k to k + 1, for k > 1. 
At step k, we have decomposed //(£,•) as a sum 

m*, o= E E ^; A2 - Ai -^ •) + E ^-Ht, •), 

l</<fc AOA23...DA; AiDA 2 D...DA fc 

where the sums run over the strictly decreasing sequences of primitive submodules of (Z d )* 
(of lengths I < k in the first term, of length /c in the second term). We have 

AiAo.-.A/^i /. /-\ / ~ AiAo...A/„i / . j- 7 7 \~1 

/^ A; (*,^s)= / A a /^a ; (*, a;, ^, a7?i, . . . , rfr?z) | T "xfl Al , 

Ji?^(Ox...xfi A ;- 1 (Ox(A i > 
J<i«)x...x^- 1 (Ox§(A fe ) 

The distributions /i^^ 2 '"^ -1 have the following properties : 

(1) ^ A 2- A '-i e C ( T * T d x §( Al ) x ... x S(A/_i) x (At))) and all its x-Fourier 
modes are in A;; with respect to £ it is supported in 7^; 

(2) for every t el, jl^ 2 "' ■) is invariant under the flows 0^ (j = 0, 1, ...,/ — 1) 
defined by 

$ £, ^) = (ar + sdH(£),£, 771, 77j_i, r?j); 

(H{ x ,^Vi,-,Vl) = + sd 2 H(^)-^-,^,ri 1 ,...,rn); 

\Vj\ 

(3) if <C then for every t Gl, /i A ; lA2 '" A!_1 (t, •) is a positive measure and 

where, for (x,£,r/i, ..,r#) G T*T d x S(Ai) x . . . x S(Aj_i) x (A,) we define: 
: f , »7i, m) 1 — ► ( x + sd 2 H(i)r] h f , r] u ...,r] l ); 
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(4) if r h = l/h then / {Aj) (t, -, d Vl ) is in C(R, M + (T*T d x§(A 1 ) x . . . xS<A,_i») 

and the measure / (Rd) * xS<Al)x ... xS<Aj _ l)x(A!) /^""^"'(t, *7i, ■ ■ ■ , dm) is an ab- 
solutely continuous measure on T d . In fact, for a £ C c °° (T*T d ) has only Fourier 
modes in A;, it admits the expression 



/ a(x, C)At lA2 "' Ai 1 ^> drji,..., drji) = 

^T*T d xS(Ai>x...xS(A i _ 1 )x(A i ) 

Tr f / aCTe i|d^ W ^.D H ~A 1 A 2 ...A ! _ 1 ( ^ ^ ; e -i^d 2 H(a)DyD y 

\ i/A 1 xS(Ai)x...xS{A ; _ 1 > 



where is the multiplication operator defined in (139]) . 2 ' 1 1 is a positive op- 
erator valued measure on J Al x S(Ai) x . . . x S(Aj_i) taking values in £}{L 2 {T d ) Aj)). 

On the other hand ^ A i A 2--- A fc satisfy: 

(1) /i A i A 2- A fe is in L°°(M, M + (T*T d x S(Ai) x . . . x S(A fc ») and all its x-Fourier modes 
are in A&; 

(2) ^ A i A 2--- A fe is invariant by the k + 1 flows, 0° : (x, £, 77) 1 — ?- (a: + sdH 771, . . . , 
and : (ac, £, 7ft, . . . , rj k ) 1 — > (x + sd 2 H(a(£))-^,£, 7ft, ... , %) (where I — 1, ... ,k). 

Finally, we define the space iS Afe which is the class of smooth functions a(x,£,?ft, . . . , me) 
on T*T d x (Ai) x ... x (A fc ) that are 

(i) smooth and compactly supported in (x,£) £ T*T d ; 

(ii) homogeneous of degree at infinity in each variable 7ft, . . . , 77^; 

(iii) such that their non-vanishing x-Fourier coefficients correspond to frequencies in A&. 

3.3. From step k to step k + 1 (k > 1). After step fc, we leave untouched the term 
£i<z< fc EaoA 23 ..oA, /ii^ 2 '"^" 1 and decompose further £ AoA23 .. oAfe /i AlA2 - Afc . Using the 
positivity of ^AiA 2 ...a^ we US e the procedure described in Section [2TT1 to write 



(40) /2 AlA2 - A *((v) = ^ ^Ai A 2-.A fc ] A 



A fc +iCA fc fc+1 



where the sum runs over all primitive submodules A^+i of Afc. Moreover, by Proposition ^. 1[ 
all the x-Fourier modes of /i AlA2 '" Afc ]^ gi? A fc ^ are in Afe + x- To generalize the analysis of 

Section 12.21 we consider test functions in <Sa +1 . We let 

c*. {.*••■■ =- (i-x(^)) 
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and 



w 



AiA 2 ...At /. _ <- \ ( Vk+l 



Rk+l 

By the Calderon-Vaillancourt theorem, both wj^*^^ Rk and tl^ 1 ^ " are bounded 
in L°°(R, («S^ 1 i ) / ). After possibly extracting subsequences, we can take the following lim- 
its : 

lim -.. lim lim(<t:l +1 (*).*) =: (A^M^a) , 

Rk+i — >-+oo Ri — >+oo h — >Q \ ' " ' K I 



and 



lim • • • lim lim 

Rk+l — >+oo _Ri — >+ooh — s-0 



*<«).«)=: (/C~**(*). 



Then the properties listed in the preceeding subsection are a direct generalisation of The- 
orems |23] and [275] (see also j3], Section 4) and writing 

(41) £ AlA2 "- A HUL ei A w = / t AlA ^ A ^(t,.,d Vk+1 )] , k 

nkeR Ak+i (a) J(A k+1 ) rikeR H+i {a > 



+ 1 ^ k+ r k (t,;dV k+ i)\ heR H {ay 



Remark 3.1. By construction, if A k+ i = {0}, we have £i AiA2 - a m-i = an d tfo e induction 
stops. Similarly to Remark \2.8\ one can also see that if rkA/c +1 = 1, the invariance 
properties o//2 AlA2 - Afe+1 imply that it is constant in x. 

Remark 3.2. Note that in the preceding definition of k-microlocal Wigner transform for 
k > 1, we did not use a parameter 5 tending to as we did when k = in order to isolate 
the part of the limiting measures supported above R^ k i (a) . This comes directly from the 
restrictions made in fijty and ( [^7] ]. 

3.4. Proof of Theorem 11.81 This iterative procedure allows to decompose /i along de- 
creasing sequences of submodules. In particular, when ~ l//i, it implies Theorem 11.81 
Indeed, to end the proof of Theorem ll.8[ we let 

Mv) = £ £ 4 lA2 - A *(v) 

0<k<d AiDA2D--OA fe DA 



= V] V] j . . ^i lA2 '" Ak (t,-,dr] u ...,dr] k )] TdxR 

0<fc<dA 1 DA 2 D-OA fc DA- / <(Ox...x<HOx(A> 

where Ai, . . . , A* run over the set of strictly decreasing sequences of submodules ending 
with A. We know that //^ lA2 "' Afe is supported on {£ G IaA, and since A C Ai we have 
Ja x C I a- We also let 

PaO)= Yl Yl a a Pa iA2 ' A * (a,dr] u ...,dr] k )] aeRAi , 

J^(f)x...x<*(0 



0<k<d AoA 2 D- 
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where the p A x 2 "' k are the operator-valued measures appearing in §3.21 

As already mentioned, Theorem 11.81 implies Theorem 1 1.2 1 in the case ~ 1/h. The proof 
of Theorem 11.21 in the case <C 1/h is discussed in Section H] and in the case ^> 1/h, in 
Section [5j 

4. Some examples of singular concentration 

4.1. Singular concentration for time scales <C 1/h. In this section, we focus on 
the case r/i < and prove Theorem 11.2( 1). 

Consider p G S (lR d ) with HpH^a^ = 1 and such that the Fourier transform p is 

compactly supported. Let (x ,C, ) G M. d x M. d and (e^) a sequence of positive real numbers 
that tends to zero as h — > + . Form the wave-packet: 

Define 

u h := Pv h , 

where P denotes the periodization operator Pv (x) := ^2 keZ d v (x + 2nk). Since p is rapidly 
decrasing, we have ||u/i|| L 2^ T d^ — > 1- It is not hard to check that (uh) is /i-oscillatory. 

Theorem 11.2( 1) is a consequence of our next result. 

Proposition 4.1. Let (77J be such that lim/ l _ > .o+ hr^ = 0; suppose that Eh 3> hr^. Then 
the Wigner distributions of the solutions S^Uh converge weakly-* in L°° (R; V (T*T d )) to 
£*(zo,£o)> defined by: 

(43) / a (x, /i(, ,&) (dx, d£) = Km i / a (x + tdff (&) , fo) Va G C c (T*T d ). 

Proof. Start noticing that the sequence (uh) has the unique semiclassical measure po = 
8 Xo <8>5f . Using property (4) in the appendix, we deduce that the image p of p (t, •) by the 
projection from T d x IR d onto M. d satisfies: 

Since for every primitive module AcZ d the positive measure JI\ is supported on R\, and 
these sets form a partition of M. d , we conclude that ~pX = unless A = A^ and therefore 
p = Pa 6o - Therefore, in order to characterize p it suffices to test it against symbols with 
Fourier coefficients in A^ . Let a G (T*T d ) be such a symbol; we can restrict our 
attention to the case where a is a trigonometric polynomial in x. Let (f G L 1 (R). Recall 
that, by Lemma |2.10[ the Wigner distributions Wh (t) of Uh satisfy 

/ ip (t) (w h (t) ,a}dt= / if (t) (w h (0) , a o <p Tht ) dt + o (1) ; 
Jr Js. 
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moreover the Poisson summation formula ensures that the Fourier coefficients of Uh are 
given by: 



u h (k) 



(27T) 



-i(k-i /h)-x Q 



Combining this with the explicit formula (1521) for the Wigner distribution presented in the 
appendix we get: 

(44) 



<p it) (w h (t) , a) dt 



(2vr 



3d/2 



k-jeA 



«0 



v[r h dH[h 



k + j 



(k - j) aj_ k h 



k + J 



£h 



h 



P ( ^ (hk - Co) ) p ( % (hj - Co) )e-^*-^ + o (1) 



-i(k-j)-x 



Now, since k — j E Ag we can write 



(k-j) 



dH h 



k + j 



< C 



h ^ — & 



-dH{So) 
\k~j\ ■ 



(k-j) 



By hypothesis, both p and k i — > at (£ ) are compactly supported, and hence the sum (J44 
only involves terms satisfying: 



h 



<R, £ ± 
h 



2 



< R and \j - k\ < R 



for some fixed R. This in turn implies 



2 ^ 



This shows that the limit of ( l44l) as /i — > + coincides with that of: 



(27T) 



3d/2 



fc-jeA 5o 
= 0(0) (0),a), 
which is precisely: 



£ 0(0) 0i _* ( ^ ) p (| - 6)) p (| (*j - 



0(O)a(x o ,£ o ) = ^(O) lim - 



a (x + tdH (f ) , fo) dt, 



□ 



since a has only Fourier modes in A^ . 

We next present a slight modification of the previous example in order to illustrate the 
two-microlocal nature of the elements of Ai (r). Define now, for i] G M. d : 

u h (x) = P [v h (x) e *»/(^)] , 
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where Vh was defined in (|42p . 

Proposition 4.2. Suppose that lini/ l _ >0 + ^ r /i — and Eh S> Suppose moreover that 

d 2 H (£ ) «s definite and that rj G (A^ ). T/ien t/ie Wigner distributions of S^Uh converge 
weakly-* in L°° (R; £>' (T*T d )) to &e measure: 

/i (t, ■) = ^( X0 +td?H{£, ) m fy)-, t e R, 
where /j,/ Xo £ \ zs £/ie uniform orbit measure defined in ( f^3| ). 

Proof. The same argumenta we used in the proof of Proposition 14.11 gives /i = /iA 5o • We 
claim that u>/ A ^ (0) converges to the measure: 

^A 5o (0, X, f , 7]) = /i(x ,£o) (^) • 

Assume this is the case, then Proposition 12.121 implies: 

&A (0 (t, X, f , 7?) = V( xo+t d?H(to)rio,ta) ( x , <^o 0?) , Vt € R, 

and, since J^a^ (t, •) are probability measures, it follows from Prop osit ion 12 . 3 1 1 hat /I A «o = 
and: 

Let us now prove the claim. Set 

u h {x)=v h (x) e ivo/{hTh) . 

Consider h > and % G C^°(lR d ) such that x^h — Uh for all h G (0, h ) and Py 2 = 1. 
We now take a G S\ and denote by a the smooth compactly supported function defined 
on M. d by a = x 2 a. Using the fact that the two-scale quantization admits the gain hr^ (see 
Remark 

(u h , Op^ (a)u h } L 2 {Td) = (u h , Op^ f0 (5)^)^2^) 

= (u h , Op^°(a)u h } L 2 {Rd) + 0(hr h ). 

Therefore, it is possible to lift the computation of the limit of Wi A ^ t h,R (0) to T*M. d x (A^ ) 
and, in consequence, replace sums by integrals. A direct computation gives: 



(u h , OpJ* {a)u h ) L , m = (2n)- d / e^^p(x)p(y) 



(x ~\~ y 1 h 1 h 

x + £h—z— > fo H H f , w(£o H ^7o H ) dxdyd£. 

Note that if F(£) = (a,rj), then 

VA; G A, F(f + fc) = (a, + fc) = F(£) + (0, k), 
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which implies that dF(C,)k = (0, k) and drj(£)k = k for all k G A^ . We deduce d7](C,Q)i] = rjo 
since £ (A&>) an< ^' ^ n y i ew of ^(£0) = 0, a Taylor expansion of r](C,) around £ gives 

ThV ( £o + — Vo + — U = Vo + o(l). 



Therefore, as /i goes to 0, 



A 



(Mh,Op h °(a)u h ) -»■ a(a;o,^o,?7o) = (j«A fo ,a). 

□ 

4.2. Singular concentration for Hamiltonians with critical points. We next show 
by a quasimode construction that for Hamiltonians having a degenerate critical point (of 
order k > 2) and for time scales <C l/7i fc_1 , the set .A4 (r) always contains singular 
measures. 

Suppose £ £ K d is such that: 

dH (f ) , (f ) , d fe-1 i7 (£„) vanish identically. 

k 

The Hamiltonian -f/"(£) = |£| (A; an even integer - corresponding to the operator (— A)a) 
provides such an example (with £ = 0). Let Uh = Pi>h, where Vh is defined in f)42p . If 
6h^> h it is not hard to see that 

\\H (hD x ) u h -H (£ ) u h \\ L2{jd) = O (h k / (e h ) k ^ . 

Therefore, 

\s{u h - ^ mo) n h \\ = tO (/i*- 1 / (e,) fc ) , 
and, it follows that, for compactly supported <p G L 1 (IR) and a G (T*T d ), 
ip(t)(w h (t) , a)dt = / ip(t) (u h , Op h {a)u h ) L2{Td) dt + O [T h h k ~ 1 / (e h ] 

Choosing (e^) tending to zero and such that ^ ^> (t/,/^ -1 ) 1 k , the latter quantity converges 
to a(xo, ^OIMIl^r) as h — > + . In other words, 

dt(g)5 X0 <gx%, G M(t), 

whence dt ® G .M(t). 

Remark 4.3. In the special case of H(£) = |£| fc (k an even integer), we know that the 
threshold is precisely h 1 . From the discussion of $3] and previously known results 
about eigenfunctions of the laplacian, we know that the elements of M. (r) are absolutely 
continuous for Th ^> . In the case of = \ jh k ~ x , one can still show that elements of 

M.(t) are absolutely continuous. This requires some extra work which consists in checking 
that all our proofs still work in this case for = 1/h 1 and £ in a neighbourhood of 
£o = ; replacing the Hessian d 2 H(!; ) by d k H(^ Q ), and the assumption that the Hessian is 
definite by the remark that [<i fc iJ(£ )-£ fc = ==>■ £ = 0] . 
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In the general case of a Hamiltonian having a degenerate critical point, the existence of 
such a threshold, and its explicit determination, is by no means obvious. 

4.3. The effect of the presence of a subprincipal symbol. Here we present some 
remarks concerning how the preceding results may change when the Hamiltonian H(hD x ) 
is perturbed by a small potential h^V(t,x). Suppose V G L°° fR x T d ) and define 

P P)h := H(hD x ) + h p V(t,x), with > 0, 
and denote by Sp h the corresponding propagator (starting at t = 0): 

Let us fix a time scale r = (77J that tends to infinity as h — > + . Define AApy ( r ) to be 
the set of accumulation points of the time-scaled Wigner distributions 

°/3,h U h 

as (uh) varies among all normalised sequences in L 2 (T d ). For the sake of simplicity, from 
now on we shall fix the time scale = 1/h. The discussion that follows can be easily 
adapted to more general time scales by changing the ranges of values of 0. 

(1) (3 > 2. In this case it can be easily shown that and Wh have the same weak-* 
accumulation points in L°° (R; V (T*T d )). Therefore, the potential is a negligible 
perturbation and, in particular, for every V G L°° (R x T d ) , 

Mp, v (l/h)=M(l/h). 

(2) (3 = 2. When H(£) = |£| 2 , the question has been addressed in [3jE3 It turns out 
that whenever V is not constant, 

M 2y (l/h)^M(l/h). 

In fact, the structure of M.2,v (VM is similar to that of Ai(l/h), but the prop- 
agation law that replaces f lTU]) involves the propagator associated to the averaged 
Hamiltonian |£| 2 + (V) A (i, x). 

(3) j3 < 2. In this case, it is possible to find potentials V for which Theorem II. 2( 2) fails, 
i.e. such that there exist \x G M.$y (1/h) such that the projection of \x on x is not 
absolutely continuous with respect to dtdx. The following example is due to Jared 
Wunsch. On the 2-dimensional torus, take H (£) = |£| 2 and V(x\,X2) '■= W(x2) 
such that W(x 2 ) = (x 2 f /2 in {\x 2 \ < 1/2}. Take s G (0, 1) and 

where \ is a smooth function that is equal to one in \\x 2 \ < 1/4} and identically 
equal to in {l^l > 1/2}- One checks that 

(-h 2 A + h 2{1 - £) V -l)u h = h 2 - £ u h + 0(h°°). 
n In that work, it is assumed that the set of discontinuity points of V has measure zero. 
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It follows that for ip G L 1 (R) and a G C c °° (T*T 2 ), 

dt 



j R ^ ( S 2(le),h U ^ °^ a ) S t-^ Uh ) L2(T2) 



= lim / ip(t) (u h} Op h (a)u h } L 2 (T 2)dt = [ / tp(t)dt] / a (x, £) fi (dx, d£) , 
h ^ 0+ Jr v ' \Jr J Jt*t 2 

and it is not hard to see that // is concentrated on {X2 = 0, £1 = 1,^2 — 0}. In 
particular the image of \i by the projection to T 2 is supported on {x2 = 0}. 

5. Hierarchies of time scales 

The following results makes explicit the relation between the sets Ai (r) as the time 
scale (r/j) varies. 

Proposition 5.1. Let (r/J and (ah) be time scales tending to infinity as h — > + such 
that lim^_ >0 + cr /i/ T "/i = 0. Then for every \i G M. (r) and almost every t 6 I t/iere exzsi 
/j* G Conv .M (cr) snc/i £/ia£ 

(45) /x (t,-) = / fJ t (s,-)ds. 

Jo 

Before presenting the proof of this result, we shall need two auxiliary lemmas. 
Lemma 5.2. Let (ah) be a time scale tending to infinity as h — > + . Let (v^ 



,(») 

y lu wbjvitwy u,s it r u . uci I C 

be a normalised family in L 2 (T d ) and define: 

T h* ,(«)■ 



h>0,nGN 



w { h n) (V) :=w[ 



Let c h > 0, n G N, be such that X]ngN c i™^ = ^-Then, every weak-* accumulation point in 
L°° (R; V (T*T d )) of 

(46) Yl c ^: ] (*. •) 



belongs to Conv .M (cr) . 

Proof. Suppose (146 j) possesses an accumulation point jl G L 00 (M; A4+ (TT*)) that does 

not belong to Conv M. (a). By the Hahn-Banach theorem applied to the convex sets {jl} 

and Conv Ai (a) we can ensure the existence of e > 0, a G C c °° (T*T d ) and # G L 1 (R) such 
that: 



and 



(£) (/i (t, ■) ,a}dt < -s < 0, 

(47) 1 6{t){fj,{t,-),a)dt>~, VfiEConvM(a). 

Jr 3 
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Suppose that fi is attained through a sequence (hk) tending to zero. For k > k big enough, 

[0(t)J2^(w^(t,-),a)dt<- 3 -e, 

which implies that there exists n fc G N such that: 

(48) / 9(t)(w^\t,-),a)dt 

Jr n 1 

Therefore, every accumulation point of (w^j also satisfies (I48jl which contradicts (I47jl . 



3 



□ 

Lemma 5.3. Let r, o and fi be as in Proposition I5.il For every a < (3 there exists 



fi>a,p Conv Ai (cr) such that 

1 



(3 — a 



H(t,-)dt = I n a% p{t,-)dt. 



Proof. Let /i G M. (r). Then there exist an /i-oscillating, normalised sequence (uh) such 
that, for every 9 G L 1 (R) and every a G C c °° (T*T d ): 

lim / 9(t)(S? t u h ,Op h (a)S? t u h )dt= [ e(jt){fi{t,-),a)dt 
h ^ 0+ Jr Jr 

Write Nh '■— Th/crh', by hypothesis Nh — > oo as h — > + . Let a < (3, define L := (3 — a 
and put: 

where L-^-^/iJ is the integer part of LA^. Then, 
I rP i rPN h 

- J (Sl ht u h , Op h (a)Sl ht u h ) L2{Jd) dt = j-jy- J ^ (S^ ht u h , Op h {a)S^ ht u h ) L . 2(Jd) dt 

lLN h \ , t h 



W E /" (S a h ht u h ,0 Vh (a)S^u h ) L2(Jd) dt 

h n =l , '*n-l 

i L^J .* h 



where the functions i>j^ := S^ htn Uy l form, for each n G Z, a normalised sequence indexed 
by /i > 0. The result then follows by Lemma 15.21 and using the fact that Sh — > 1 as 
h — >0 + . □ 

Proof of Proposition 15. 11 Let fi G Ai(r); an application of the Lebesgue differentiation 
theorem gives the existence of a countable dense set S C (T*T d ) and a set N C R of 
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measure zero such that, for a G S and t G R \ N, 

(49) lim — / / a (x, £) fi (s,dx,d£) ds = f a (x, £) \i (t, dx, d£) . 

e^0+ 2e J t _ £ J T *T d JT*T d 

Fix t G R \ N; then, for any e > there exist /i* G Conv (a) such that, for every 
a G C c °° (T*T d ) , 

(50) — — / / a(x,£) [i(s,dx,d£)ds = I I a (x, £) //* (s, dx, d£) ds. 

2 £ Jt-e JT*T d JO JT*T d 

Note that ConvA / l(o") is sequentially compact for the weak-* topology, therefore, there 
exist a sequence (e n ) tending to zero and a // G Conv .M (cr) such that /i* converges 
weakly-* to //. Identities (|4"9]) and (150)) ensure that yU (t, ■) = J (s, ■) ds. □ 

Remark 5.4. Projecting on x in identity we deduce that given v G Ai (r) i/iere exzsi 
(/' 6M (a) such that: 

v (t, •) = / z/' (s, •) ds. 



o 



77ms 7 together with the fact that elements of Ai (1/h) are absolutely continuous imply the 
conclusion of Theorem M .2( 2) when ^> 1/h. 

Denote by Ai (oo) the set of weak-* limit points of sequences of Wigner distributions 
{w Uh ) corresponding to sequences (uh) consisting of normalised eigenfunctions of H (hD x ). 
We now focus on a family of time scales r for which the structure of Ai (r) can be described 
in terms of the closed convex hull of Ai {oo). Given a measurable subset O C R d , we define: 

T/ f (O) := hsup{\H (hk) - H (hj)\~ X : H (hk) + H (hj) , hk, hj G hZ d n 0} . 

Note that the scale defined in the introduction coincides with (R d ) . The following 
holds. 

Proposition 5.5. Let O C R d be an open set such that (O) tends to infinity as h — > 

+ . Suppose (r/j) is a time scale such that lim/ l _ s>0 + 7/f (O) /r^ = 0. If \i G Ai (r) is 
obtained trough a sequence whose semiclassical measure satisfies //„ (T d x (R d \ O)) =0 

i/ien G Conv.M (oo). 

Proo/. As in (PJ, for a G C~ (T*T d ) and 6 £ L 1 (R), we write: 
(t) (iu A (t) ,a)dt = 

1 7T H (hk) - H (hj)\ ^ A,, .A 

1 ; hjezd v 7 v 7 
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Our assumptions on the semiclassical measure of the initial data implies that, for a.e. 
t E R: 

fi(t,T d x (R d \0)) =0. 

Suppose that fi is obtained through the normalised sequence (uh)- Suppose that a G 
Cf (T d x O) and that supp# is compact. For < h < ho small enough, 

^ H(hk)-H(hj) ^ gupp ~ y hk ^ h . ^ Q guch that H ^ ^ R ^ ^ 



Therefore, for such h, a and 6, 



(t) (w h (t) , a) dt - ° { ° ' 



v ' kh,hjeO 

H(hk)=H{hj) 

= 0(0) Yl (PE h u h ,Op h (a)P Eh u h ) L2m , 

E h £H(hZ d )nH(0) 

where stands for the orthogonal projector onto the eigenspace associated to the eigen- 
value Eh- This can be rewritten as: 



/ 9(t)(w h (t),a)dt = 9(p) V 



E h £H(hW)C\H(0) 

where 

E h ._ PE h Uh , E h ._n p i|2 

\\PE h u h \\ L 2^ \ > 

Note that v^ h are eigenf unctions of {hD x ) and the fact that is normalised implies: 

£ = i- 

E h eH(hZ d )nH{0) 

We conclude by applying (a straightforward adaptation of) Lemma to v h h and c^ h . □ 

Corollary 5.6. Suppose := (R d ) — oo as h — > + and i/ia£ (r^) is a time scale 
such that <C t^. T/ien 



.M (r) = Conv A4 (oo). 



Proof. The inclusion .M (r) C Conv M. (oo) is a consequence of the previous result with 
O = M. d . The converse inclusion can be proved by reversing the steps of the proof of 
Proposition 15.51 □ 

Remark 5.7. Proposition \ 1 . 1J\ is a direct consequence of this result. 
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6. Appendix: Basic properties of Wigner distributions and semi-classical 

measures 

In this Appendix, we review basic properties of Wigner distributions and semiclassical 
measures. Recall that we have defined w\ for Uh G L 2 (T d ) as: 



(51) / a(x, W L (dx, d£) = (u h , Op h (a)u h } L2(Jd) , for all a G C™(T*T d ), 

JT*T d ' ' 



Start noticing that pTJ) admits the more explicit expression: 



1 ^ fh 



(52) / a(x,£)w Uh (dx,d£) = (2 y /2 2^ u h {k)u h {J)aj- k ( ^(k + J 

where u h (k) := f Jd u h {x) *~* d/2 dx and a k (g) := f Jd a(x,£) ** d * 2 dx denote the respective 

Fourier coefficients of Uh and a, with respect to the variable x G T d . 

By the Calderon-Vaillancourt theorem [8], the norm of Op ft (a) is uniformly bounded in ft,: 
indeed, there exists an integer K^, and a constant > (depending on the dimension d) 
such that, if a is a smooth function on T*T d , with uniformly bounded derivatives, then 

||Op 1 (a)|| i 2( T d)_ >L 2 (P t) <C d ^ SU P l^° a l =: C dM{a) . 

aeN 2d ,\a\<K d T * Jd 

A proof in the case of L 2 {R d ) can be found in [TD]. As a consequence of this, equation (j5"Tj) 
gives: 



< C d |K||' 2(Td) M (a) , for all a G C(T*T d ). 



Therefore, if Wh{t, •) := for some function ft i — > 7} t G R + and (it^) is bounded in 

one has that (wh) is uniformly bounded in L°° (R; £>' (T*T d )) . Let us consider 
/it G L°° (R; X>' (T*T d )) an accumulation point of (wh) for the weak-* topology. 

It follows from standard results on the Weyl quantization that fi enjoys the following 
properties : 

(1) \x G L°°(R; Ai+(T*T )), meaning that for almost all t, [i(t, •) is a positive measure 
onTT d . 

(2) The unitary character of implies that j T , Td fi(t, dx, d£) does not depend on t; 
from the normalization of Uh, we have j TtJd fi(r, dx, d£) < 1, the inequality coming 
from the fact that T*T d is not compact, and that there may be an escape of mass 
to infinity. Such escape does not occur if and only if (uh) is /z-oscillating, in which 
case /i G L°° (R; V (T*T d )). 

(3) If Th — > oo as h — > + then the measures •) are invariant under <f> s , for almost 
all t and all s. 

(4) Let Ji be the measure on R d image of fj,(t, •) under the projection map (x,£) i — > £. 
Then /2 does not depend on t. Moreover, if /Z^ stands for the image under the same 
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projection of any semiclassical measure corresponding to the sequence of initial data 
(u h ) then Ji = Jj^. 

For the reader's convenience, we next prove statements (3) and (4) (see also j20] for a 
proof of these results in the context of the Schrodinger flow e lhtA on a general Riemannian 
manifold). Let us begin with the invariance through the Hamiltonian flow. We set 

a s (x, £) := a(x + sdH = a o <f) s (x, £). 

The symbolic calculus for Wey's quantization implies: 

±S' h Op h (a s )S~ s = S s h Op h (d s a s )S^ s - j-SJl [H{hD) , Op h (a s )\ S h s 
= 0{h% 

Therefore, S s h Op h (a s )S^ s = Op h (a) + 0(h 2 ) and for 9 G L X (R), 



6(t)(w h (t),a)dt= I 9(t)(u h , S- Tht Op h (a)S T h ht u h )dt 

B(t)(u h , S- Th(t - S/Th) Op h (a o 4> s )S T h h{t - s/Th) u h )dt + 0{h 2 ) 
6(t + s/r h )(u h , S~ Tht Op h (a o cj> s )S Tht u h )dt + 0(h 2 ) 
6(t + s/r h ) (w h (t),ao <p s )dt + 0(h 2 ). 

Since \\9{- + s/rh) — 0\\li — > (recall that we have assumed that r^ — > oo as h — > H 
we obtain 

6{t) (w h (t) ,a)dt- / 0{t)(w h (t),ao <j> s )dt — > 0, as h — > + , 



whence the invariance under 

Let us now prove property (4). Consider ~fi the image of /i by the projection (x, £) i — > £, 

IOC 

-0 



we have for a G C?? 1 ™^ 



-t d 

K(t),a(0)-«,a(0>= / -r-iWhis) , a{C))ds 



ds 

1 d 
Uh , -J- (S h ThS Op h (a)S T h hS ) u h )ds 



o 

= 0, 

as 4^S^ S Op ft (a(^))S'^ = (for a only depending on £ we have Op h (a) = a(hD x ), which 
commutes with H(hD x )). Therefore, taking limits we find, for every 9 G L 1 (R): 

6 (t) I a (0 pl (t, dx, <%)=( fe (t) dt] ! a (0 no (dx, d£) , 

JT*T d \JR J JT*T d 
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where /i is any accumulation point of (w^J . As a consequence of this, we find that \x 
does not depend on t and: 

7*(0 = / Vo{dy,Q. 
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